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We argue that the gauge invariant source terms in the formalism of the Color Glass Condensate 
(CGC) should properly describe the target wave-function if the microscopic structure is concerned in 
the regime where the parton distribution is not quite dense. The quantum property of color charge 
density is incorporated in the quantum weight function defined with the source terms. We sketch 
that the logarithmic source terms encompass a meaningful picture of the microscopic structure of 
the target wave-function. 
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\£) • Introduction The formalism of the Color Glass Condensate (CGC) is elaborated suitably to describe energetic 
matter in Quantum Chromodynamics (QCD) 0,0, When the target material like a nucleus whose dynamics is 
mainly governed by QCD is highly boosted along the positive z direction, more and more partons with small fraction 
of the longitudinal momentum p + — (p° + p z )/-i/2 take part in the dynamics. Such partons are commonly called the 
wee partons. The basic idea of the CGC is as follows; the wee partons with p + smaller than a certain separation scale 
I ^ ' A distribute densely at sufficiently high energy that the wee parton distribution is well coherent and can be described 
by the classical field A%(x). The description is analogous to classical electromagnetism then once the saturation is 
reached. The wee parton distribution is to be approximated by solving the classical equations of motion, i.e. Yang- 
Mills-Maxwell equations. The color source p a (x) which generates A£(x) is provided by fast partons, namely, partons 
£SJ ■ with p + larger than A. 

\ The name of the CGC is widely understood: Color is the quantum number relevant to QCD partons. Glass is not 
crystalline microscopically but static enough as compared with our daily time scale. Likewise, the fast-moving partons 
are relatively static seen from the wee partons due to the time dilation. Condensate means that the wee partons are 
dense and coherent represented by a classical field like the field-theoretical description of the Bose-Einstein Condensate 
(BEC). We note, however, that the physics of the CGC is totally different from that of the BEC. The CGC is not a 
, state of matter, but rather an effective description of matter. 

The separation scale A is physically specified by Bjorken-x which intuitively signifies the longitudinal momentum 
fraction between the target's P + and the parton's p + in the infinite momentum frame where P + ~ oo. The random 
distribution of p a (x) depends on A or x through the weight function Wa [p] . The differential equation of Wa [p] with 
respect to A is the evolution equation by which one can go toward smaller A and thus smaller x. The evolution 
equation is an essential tool to investigate the small-x physics. In the leading-logarithmic approximation the CGC 
formalism leads us to the evolution equation known as the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner 
(JIMWLK) equation HSH 0, II H E3 • Although the JIMWLK equation is successful when the color charge density 
p a (x) is as large as ~ 1/g in the saturation region of Q 2 and x, it has turned out that an important process in lower 
density regimes is missing in the JIMWLK dynamics. That is the fluctuation of the parton (gluon) number in the 
target wave-function. The JIMWLK equation correctly grasps the saturation mechanism by the mergi ng p rocess of 
gluons, while it fails in describing the splitting process that is indispensable in the dilute regime [TH fl3 | . 

It is obvious what lacks in the JIMWLK dynamics. If the target material of the charge density p a (x) travels 
through the gauge field A% (x) created by the projectile, then the scattering amplitude earns an eikonal phase ~ 
exp[— ip a (x)A~ (x)], that means there is a sort of uncertainty principle between the target and the projectile states. 
Thus, if the eigenstate of the target wave- function is chosen represented in terms of p a (x), the projectile wave-function 
is uncertain represented in terms of i8/6p a (x), which is analogous to x\x) = x\x) and p\x) — ihd/dx\x) in Quantum 
Mechanics. The JIMWLK equation contains the full order in p a (x) but only the first and second order in i8/5p a (x) 
under the assumption that p a (x) is large (and thus iS/Sp a (x) is small). 

The higher order term in id/5p a (x) has been added in a way consistent with the color dipole model (DM) in 
Refs. [j^ll&IIij and the general evolution kernel or the Hamiltonian including full order in both p a (x) and i5/Sp a (x) 
has been achieved in Refs. 0, 0, 0, 0] . The state evolving according to the Hamiltonian should be specified to 
write down an evolution equation. Thus, Refs. 0,0, H(j are seminal works in which they showed that the DM weight 
function evolves properly according to the general Hamiltonian in the dilute regime (Bremsstrahlung Hamiltonian) in 
the large N c limit. 

It is not quite clear, however, how we can push the conventional strategy to derive the JIMWLK equation toward a 
lower density region. Apparently, for one example, the derivation of the general Hamiltonian in Rcf. |l7fl does not rely 
on any information of the target nor projectile wave-function, while the JIMWLK equation deals with the evolution 
of the target wave-function. Its clear manifestation is; when one derives the JIMWLK equation, one needs to have 
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the gauge invariant source terms that induce complicated and non-local vertices. In contrast, the general Hamiltonian 
shows itself without complicated vertices associated with the source action at all. This seemingly mystical simplicity 
as compared with the derivation of the JIMWLK equation is to be understood as follows 0, |23 ■ The source terms 
should carry the information of the target wave-function. The general Hamiltonian works when it acts either on 
the target wave-function that is embodied as the weight function with the source terms as we will see, or on the 
projectile wave-function that is supplied as an operator to be evaluated with the weight function. Even though the 
Hamiltonian is reduced to the JIMWLK kernel in the dense case, it is not the JIMWLK equation until the weight 
function is specified. Since the parton distribution is highly dense in the JIMWLK problem, the concrete form of 
the source terms is not relevant to the final result. In the dilute regime, however, the explicit form is significant as 
demonstrated in the choice of the weight function as in the DM in Ref. J2(j . This point of view is already implied in a 
slightly different formulation of the high-energy evolution in Refs. 0, Il9[ . while we will address here how we should 
understand this along with rather conventional JIMWLK formulations 0, 0, 0J • 

We will start with the generating functional and briefly review how the CGC picture arises from the stationary-point 
approximation. Then, we illustrate a simple algebraic trick proposed by the present author in Refs. 0, to treat 
the quantum nature of color charge density. The benefit of this trick is also that the physical meaning of the source 
terms is transparent. Actually it becomes clear that the source terms should describe the target wave-function, that 
is, in the dilute regime the source terms should comply with the microscopic structure of the target wave-function. 
We will claim that the logarithmic form would be a likely candidate. 

Stationary-point approximation In the presence of the color source that distributes with the weight function 
W r [p], the QCD generating functional is modeled as 



Z\j} = JvpW T [p] J VA5[A+] expjzS YM L4] 



+ iS w [p]-i d 4 x 3 l{x)A^x)\ (1) 



in the light-cone gauge. Here £Vm[^4] is the Yang-Mills action and Sw[p] denotes the source action which is given in 
terms of the gauge invariant variable, 



W[A~](x~,x) = V x + exp|i 5 y dx + A~ (x + , x~ ,x) j , (2) 

when the background field comes only from the target moving along the x + direction. It should be noted that A~ 
in the above expression is a color matrix either in the fundamental or adjoint representation in accord with the 
constituent of the target. It is the central issue of our interest how we should fix the concrete form of SV[A - ]. In 
writing Eq. we introduced the rapidity variable, r = ln[P + /A] = ln[l/x], instead of the separation scale A or 
Bjorken-x. Then, one can compute the expectation value of an operator 0[A] as 

{0[A]) T = 0[i8/5j]Z\j]\. = jZ[Q]. (3) 

The CGC results from the stationary-point approximation on Eq. (JJJ. When p a (x) is as large as ~ 1/g at small 
x, the classical approximation works well, so that the functional integral over A^ can be estimated at the stationary 
point determined by 

Then the expectation value ® is approximated as 

(0[A]) r = JvpW T [p]0[A[p\] (5) 

under the assumption that 0[A] would hardly affect the stationary condition. In case of computing the scattering 
amplitude between a heavy target of the CGC and a simple projectile like the color dipole, the above stationary-point 
approximation is acceptable, while one needs to solve the two-source problem at the classical level if one intends to 
compute observables like the gluon or quark production associated with the scattering [23, l2 ll l25l l2q | even though it 
is a dense-dilute system. Here we only care about the target wave-function and pay little attention to the projectile. 

Let us consider what is required for the property of Sw [p] . The leading term of Sw [p] should reproduce the eikonal 
coupling — Jd 4 x p a (x)A~(x). From the equations of motion @ the charge density inferred from the source terms 
should satisfy the covariant conservation, which is fulfilled whenever the source terms are written in terms of the 
matrix elements of Eq. J5J) as one can easily confirm 

{d-6 ab + gf acb A-)j+ = 0, 

.]■ x (U\ ,.. igT b W x+> x> )... 
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Here T° is an element of the SU(iV c ) algebra. The covariant conservation holds for both the fundamental and the 
adjoint Wilson lines. In Eq. JSJ we denoted the Wilson lines with the x + integration ranging from x + to oo and from 
— oo to x + as Woo x+ an d VV x + _ oa respectively. 

In the dense regime there is no further restraint onto the choice of SV[p] and, though it is not rigorously proved 
yet, there is likely a sort of universality and the concrete form of Siy[p] is n °t relevant to the evolution equation, i.e., 
the JIMWLK equation. 

A further constraint appears in the dilute regime where the color charge distribution in the target is not quite 
dense. For that purpose it is convenient to transform the generating functional into the following representation; 



Z\j] = J VpN T [ P ] J VA5[A+] expjiSW^] - i J d 4 xp a { 



(x)A-(x)-i d*xj a Jx)A»(x) , (7) 



N T [p] = f T>pW T [p} exp{iS w [-i S/6p, p}\ S[p}. (8) 



where 



This is just a simple rewriting, and nevertheless, it can bridge the formalism toward a DM-like picture as we will see. 
Later we will closely discuss the physical meaning of the expression (jHJ. 

It is a notable feature of this representation that one can see the non-commutative nature of color charge density 
manifestly. The stationary-point of Eq. Q is determined by the equations of motion; 

SSym s ^ pa \ = Q> (g) 



and the expectation value of 0[A] is then approximated as 

(0[A]) T = fvpM T [p}0[A[p}]. (10) 



It is important to note that p a (x) is dependent on x + , which actually takes care of the full account of SSw/^A^ in 
Eq. (@J. Therefore, we can no longer assume that the color source p a {x) is static once the non-commutativity plays a 
role. In the functional integral formalism we have the commutation relation as follows; 

(p a (x+,x)p b (x+- V ,y) - p\x+,y)p a (x+- r/,f)).^ 0+ 



= VpM T [p}{p a (x + } x)p b (x + - V ,y)-p b (x+,y)p a (x + ~r 1 ,x)} 



7)->0+ 



= - (igf jVpNM ^ (W^ x+ [T*,I*\W B[ + t _ 00 )..6<n(S-$ 



sw l3 

= - igf abc ■ i (ig) J VpM T [p\ ^ (W x>x+ T c W x+ ,_ x ).. S (3 \x - y) 

= -igf abc {p c (x + ))5^{x-y). (11) 
Here we utilized the following relation, 

dxf(-i-^jS(x)x = Jdx{f(Q)-if'(0)-^----}s(x)x = if'(()). (12) 

Strictly speaking, in fact, the right hand side of Eq. I|ll[) is zero after average. One can insert an appropriate operator 
to make the expectation value nonvanishing and thus Eq. I|ll|l should be understood in such a sense. What we have 
seen here is essentially equivalent with the argument by Kovner and Lublinsky in Ref . [34| , though they just mentioned 
a similar form encountered in the DM as a motivation. We emphasize that the argument so far is applicable for any 
SV[/o] as long as SV[p] is a function of the Wilson line in the previous works 0,H3 the present author called 
M T [p] the density of states. Here it would make sense to call Af T [p] the quantum weight function since it properly 
handles the quantization relation. 

Let us shortly comment on a symmetric representation of the scattering problem of light materials. In that case, 
as the target is a right mover, we shall denote the above p a {x), W T [p], and <SV[p] as Pr( x )> VV^'[/0r], and S'^/Ipr] 
specifically. In general the wave-function of the left-moving projectile can be expressed as 

0[A] = J Vp L W£_ T [p L ] exp{^L4+ p L ]} (13) 
PpL^PL Wy-M exp(i£^HV<JpL,PL] - i ! d 4 xpl(x)A+(x)\ 5[p L ] (14) 
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with Y being the total relative rapidity of the target and projectile. Then, if the scattering amplitude is concerned 
and A + in Eq. II Kit! is simply approximated as -4 + [p] like in Eq. I|1(J|I . the scattering amplitude is evaluated as 

S Y = J V P Vp K Vp h W T R [p R ] W y _ t [pl] ex.p{iS^[-i6/5p, p R ] + iS^[A\p], pl}} 5[p}. (15) 

This expression makes sense in the dilute-dilute or dense-dilute (where the dense-dilute duality |27| realizes) scat- 
terings, but not in the dense-dense case in which the classical solution A^ [p] must be evaluated in the presence of 
two dense sources of the target and projectile. The color glass representation of the onium-onium scattering like as 
formulated in Ref. [3^ | is one manifestation of the above. If one is interested in an expectation value of some other 
operators, such as in the gluon production pEj and in the quark production |26l | and so on, or if one is interested in 
the dense-dense case, one has to use Eq. (|14[) and solve the two-source problem with pn(x) and pl{x) as analyzed in 
Refs. [23, |2J, uli 1211 l2a 1211 Lttl LUJ • Now let us return to the main stream of this work and discuss the choice of the 
source terms. 

Choice of the source terms We know that in the large N c limit the DM weight function captures the correct 
description of the target wave-function in the dilute regime. The dressed dipole weight function reads 

00 . N 

Z T DM [p] = W dT N PB M ({z t }-r)l[DHz^ 1 ,zi)S[p}, (16) 

JV=1 J i=l 

where N is the number of dipoles and dTjy means the integration over all the phase space, namely, dT^ — dz\- ■ ■ dzjss 
and Z)t (x, y) is the creation operator of a dressed dipole located at x and y in transverse space, that is, 

D^{x,y) = ^tT[W{x)W\ V )] (17) 

in the fundamental representation. In the DM the longitudinal structure in the x~ direction is neglected and the 
Wilson lines do not depend on the longitudinal coordinate. Now one may well notice that Eq. (|16fl is a special case of 
the quantum weight function JSJ. In other words, the source terms SV^ - ] must be chosen such that the quantum 
weight function can reproduce Eq. (| 1 fc>|) somehow in the large N c limit; 

N T [p] — ► Z? M [p] as iV c ^oo. (18) 

It is important to mention that we always mean the dressed dipoles written in terms of the Wilson line. Actually, in 
the weak field limit, any source terms which satisfy SW[-<4] = — J d 3 x p a (x)A~{x) + ■ ■ ■ would lead to the results of 
the (undressed) DM if only an ensemble of the dipole-type color source, 

p a (x) = Q a 8(x-){S^(x - B ) ~ 5^ ( X - „„)}, (19) 

where Q a is a Gaussian random variable ^3 > is allowed by W^foj. 
The conventional choice adopted in literatures 0, 0, @, S El EJ 



Sw[p] = ^- J d 3 xtv{p(x) W[A-](x)}, (20) 

does not meet the requirement of Eq. (|18|) : whatever W T [p] cannot lead to the form of the weight function (flfjfl . 
Therefore, the conventional form of the source terms is no longer a right choice microscopically and is ruled out. It is 
instructive to look further into why the conventional choice is not good. If the above conventional source terms are 
substituted into the weight function © and the exponential is expanded, an arbitrary number of W is allowed at each 
point. The higher order W terms cannot be controlled by the distribution W T [p] and they are not compatible to the 
specific form of l|ltj|) . To put it in a more formal way, this incompatibility seems to stem from that the color charge 
density belonging to the algebra and the Wilson line belonging to the group are put on the equal footing, which is 
quite unnatural. 

One possibility would be the eikonal coupling with the Wess-Zumino action [3^1 Isfij that takes care of the 
non-commutative nature of color charge density just like our M T [p\. The Wess-Zumino action Swz[p] is a function of 
p a {x). In principle 5wz[p] is given in terms of W[— iS/Sp], but in fact, W is implicit in the expression of Swz[p] given 
111 Refs. and it is non-trivial how to retrieve W out of <Swz [p] ■ We have no idea how to derive the microscopic 

information of the target wave- function, not going back to the representation as a product of W but directly using 
Swz [p] on its own. 
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Instead, we shall here focus on another specific possibility of the source terms. We will check whether the following 
logarithmic form jM El HI , 

(21) 



Sw[A-,p] = i- I d 3 xtr{p(x)lnW[A-]{x)}, 



in the fundamental representation, or 



S w [A-,p]=i— I d 3 xtx{p(x)lnW[A-](x)}, 



(22) 



in the adjoint representation (W is defined with the adjoint generators, and so is p(x) implicitly), is qualified or not 
as a microscopically meaningful choice. In this form we note that In W belongs to the algebra unlike W , so that the 
coupling to p a (x) looks reasonable. 

We shall develop one more argument besides Refs. 0, 122I Isfil ] that the logarithmic form (|21|) is motivated. We 
know that the solution of the covariant conservation (jfJJ takes a form of 



j:(x + ,x) = w^^ 00 p a _ oa 

with the initial condition being p a _ x at infinite past. Using the relations W ha — 2ti\t b Wt a W^\ and WlL 
^00 -00 ^oo,x+ 1 we can equivalently rewrite the above as 



(23) 



2 

-t — tr 
9 



P- 



sw 

5A7 



i — tr 
q ^ 



P-ooil-W) 



J (i-w) 

SAZ 



(24) 



where we expanded W' = {1 — (1 — W^)} -1 in terms of 1 — W . Because 1 — W is at least linear in A~ if expanded, 
intuitively, (1 — W) n in the above represents the (at least) n gluon insertion. Now we change the variable p-00 so as 
to allow for all the orderings of the source among the gluon insertion, namely, 



P-oo(l-W) n - 

Then, Eq. simplifies as 



-lp (1 - W) n + (1-W)p(l- WY 1 - 1 + (1 - Wf p (1 - W) n - 2 + •••}. 
n l > 



J h + = -i -tr 



SlnW 



sa; 



(25) 



(26) 



which exactly coincides with the current deduced from the logarithmic source action (|21|l . From this argument, it is 
also clear that our p(x) does not correspond to the initial condition at infinite past. We did not derive Eq. l|21fl from 
Eq. J^J, for Ea. H25|) is a highly non-trivial variable change, but we can insist that Eq. I|21() rather than the conventional 
one is well motivated. 

Returning to the main discussion, let us draw an intuitive picture of the microscopic description of the target wave- 
function which consists of partons propagating along the x + direction. The classical source p a (x) indicates the parton 
distribution. It is quite important to notice that the distribution of the color orientation of p a (x) is not random at 
each spatial point. Of course, in the dense regime, one can regard the coarse-grained color distribution as random at 
each point, but one cannot in the dilute regime. Otherwise there would be no correlation between different two points 
of a quark and an antiquark and the dipole-type operator like Eq. I|17|l is not permitted but only trW (x) trW^ (y) 
would arise. The random distribution is expected to be valid at scale longer than the saturation scale ~ Qj 1 and the 
dipole size would be smaller. The above mentioned is explicitly observed in Eq. 119|) where the random variable is 
common for the quark at Xq and the antiquark at y . 

In order to see the meaning of the logarithmic form, for the simplest example, let us assume that a quark and an 
antiquark sit at Xq and y respectively in the same way as in Eq. (|19|l . (We will consider the adjoint case later.) As 
we explained above, the spatial resolution of randomness is coarser than the separation scale of the dipole. Then, the 
source part becomes 



exp 



--tr{p(\nW(x )+\nW\y ))} 



(27) 



where we used the fundamental one (|21|) and — In W — In 14^. In general, of course, there could be a color rotation 
between p a (xo) and p a (y ). If p a (y Q )t a = S p a (xo)t a , however, such a color rotation by S can be absorbed by the 



FIG. 1: Intuitive picture of the leading contribu- 
tion in the N c counting. 



FIG. 2: An example of the generalization to the 
case with multiple adjoint particles. 



redefinition S'W"* (y a )S — > W'(y Q ). The point here is that the random average over p a is taken simultaneously for 
the quark and antiquark and Eq. (|27() is then a natural extension of Eq. I|19|) to the non-linear situation. 
Here one can make use of the Campbell-Hausdorff formula, 

\nW(x )+\nW^y ) = \n[W(x Q )W^y )]~^[\nW(x ), In W^(y )] + ■ ■ ■ . (28) 

In the large N c limit only the first term is dominant; the gauge coupling scales as g ~ l/y/Nc which is compensated 
by quark loops made up in the double line notation of gluons. The higher order commutators lead to the higher order 
terms in g and yet the color structure is reduced by the commutators, and thus they are suppressed. At the price of 
mathematical rigor we sketch the intuitive picture of this in Fig. ^ Recalling that the source terms originate from 
the QCD vertex g"tjjiA~t^ipj in case of the fundamental representation, Eq. lj2T|) can be expressed as the one-gluon 
vertex between a propagating quark and gluons as in Fig. ^ Here we adopted the double line notation for gluons at 
Xq and y . Roughly the first term of Eq. (|28|> corresponds to the contribution with the dotted line in Fig. ^ that 
makes a loop and thus is the leading order in the N c counting. 

We drop the sub- leading commutators in Eq. (|28|l . Therefore, in the presence of a quark and an antiquark and in 
the large N c limit, the quantum weight function reads 



Mr\p]~ JdpW T [p]exp ~tt{p\n[W(xo)W^(y Q )]} 



5[p}. (29) 



Now the calculation goes in the exactly converse manner as in Ref. [2l|. From the property, VV r [Sp>S^] = W r [p], one 
can diagonalize W(x )W ji (y ) to have a polynomial in terms of its eigenvalues. Using W r [S'pS' t ] = W T [p] again, one 
can conclude that the polynomial is to be expressed in terms of the trace of W '(xq)W^ (y ) (see also the discussions 
in Ref. HJ). Hence, when W T [p] allows for up to the first order of the trace, we reach 

Mr\p]~Dt(x ,y )8\p]. (30) 

What we showed here is that there is a certain W T [p] that makes the weight function © with the logarithmic source 
terms l|21f) being the DM weight function (|16fl . but there is not for the case with the conventional source terms i|2U[l . 

The generalization to the case with more quarks and antiquarks is straightforward under the assumption that the 
randomness of color orientation is clustered at the scale of dipoles. Also, the generalization to the adjoint source l|22|l . 
i.e. the case with the target composed of gluons, is easy too. The above argument works as it is with the Wilson line 
replaced by the adjoint one W. The formula, t^t b kl — \8u5jk — j]y-o~ij5ki, follows that the leading order at large iV c is 
given by all the dipole-type combinations of gluon double lines, as exemplified in Fig. [5] Thus, as an ensemble, these 
combinations eventually lead to the DM picture with a certain weig ht P° M . 

Discussions We demonstrated that the logarithmic form, at least in principle, encompasses the desirable micro- 
scopic structure of the target wave-function. It had been confirmed that the logarithmic source terms correctly lead 
to the JIMWLK equation in the dense regime, and in addition to that, it has been checked here that they are qualified 
also in the dilute regime. In fact, however, the logarithmic form of the source terms alone does not provide sufficient 
microscopic information. That means, in order to determine the dynamics in the dilute regime, the microscopic 
information of W T [p] is required; especially the randomness is only at the unit of the dipole size scale, or in other 
words, the color orientation has a strong correlation at short distance of the dipole size. The coarse-grained picture 
of the color average is valid only in the dense regime. 

Let us assume that we had a complete evolution equation applicable for the dense, dilute, and intermediate regimes. 
It would be an intriguing question what would transpire at intermediate density. When the source is large enough in 
the dense case, the color distribution is regarded conventionally as random. That is, in two-dimensional transverse 
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space, it is completely disordered. When the density gets lower and lower, the finer structure of randomness becomes 
significant, and the two-dimensional system is no longer disordered and a sort of localization is important. We have 
no idea so far what dynamics would separate these two distinct regimes. Is it a continuous crossover from one to the 
other, or does it result in a discontinuous transition between two regimes? 

It is known that neither of the JIMWLK and the DM equations would be a good starting point to tackle this 
problem in the intermediate density region. The JIMWLK dynamics is valid after coarse graining and thus it contains 
less information. The DM might be good, but it is unlikely that the DM is capable of taking care of saturation, for 
which the interaction among partons is important. We would expect that our logarithmic expression may give a hint 
for future developments. 
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